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Abstract. We discuss a novel cold dark matter candidate which is formed from 
the ordinary quarks during the QCD phase transition when the axion domain wall 
undergoes an unchecked collapse due to the tension in the wall. If a large number of 
quarks is trapped inside the bulk of a closed axion domain wall, the collapse stops due 
to the internal Fermi pressure. In this case the system in the bulk, may reach the 
critical density when it undergoes a phase transition to a color superconducting phase 
with the ground state being the quark condensate, similar to the Cooper pairs in BCS 
theory. If this happens, the new state of matter representing the diquark condensate 
with a large baryon number B ~ 10 32 becomes a stable soliton-like configuration. 
Consequently, it may serve as a novel cold dark matter candidate. 
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1. Introduction 

The presence of large amounts of non-luminous components in the Universe has been 
known for a long time. In spite of the recent advances in the field ( see e.g. recent 
summary pQ), the mystery of the dark matter/energy remains: we still do not know 
what is it. The main goal of this work is to argue that the dark matter could be nothing 
but well-known quarks which however are not in the "normal" hadronic phase, but 
rather in some "exotic" , the so-called color superconducting (CS) phase. 

This is a novel phase in QCD when light quarks form the condensate in diquark 
channels, and it is analogous to Cooper pairs of electrons in ordinary superconductors 
described by BCS theory. There existence of CS phase in QCD represents our first 
crucial element for our scenario to work. The study of CS phase received a lot of 
attention last few years, see original papers j2j,[Sl and recent reviews jlj on the subject. 
It turns out that CS phase is realized when quarks are squeezed to the density which is 
few times nuclear density. It has been known that this regime may be realized in nature 
in neutron stars interiors and in the violent events associated with collapse of massive 
stars or collisions of neutron stars, so it is important for astrophysics. The goal of this 
work is to argue that such conditions may occur in early universe during the QCD phase 
transition. Therefore, it might be important for cosmology as well. 

The force which squeezes quarks in neutron stars is gravity; the force which does 
a similar job in early universe during the QCD phase transition is a violent collapse 
of a bubble formed from the axion domain wall. If number of quarks trapped inside 
of the bubble (in the bulk) is sufficiently large, the collapse stops due to the internal 
Fermi pressure. In this case the system in the bulk may reach the critical density when 
it undergoes a phase transition to CS phase with the ground state being the diquark 
condensate. These configurations with large number of quarks in color superconducting 
phase, will be named the QCD balls. Therefore, an existence of the axion domain wall 
represents our second crucial element for our scenario to work. We should note at this 
point that the axion field was introduced into the theory to explain the lack of CP 
violation in the strong interactions. Later on the axion field became one of the favorite 
candidates for the cold dark matter, see original papers [Sj-jHl and recent reviews |Hj on 
the subject. In the present scenario the axion field plays the role of squeezer rather than 
dark matter itself. In principle, it can be replaced by some other, yet unknown fields 
with similar properties. However, to be more concrete in estimates below we shall use 
the specific properties of the axion field with known constraints on its coupling constant. 

We do not address the problem of formation of QCD-ball in this letter. Instead 
we concentrate on the problem of stability of these objects. As we will show, once 
such a configuration is formed, it will be extremely stable soliton like particle. The 
source of the stability of the QCD-balls is related to the fact that its mass Mb becomes 
smaller than the mass of a collection of free separated nucleons with the same baryon 
charge. The region of the absolute stability of the QCD-balls is determined by inequality 
m N > M B — M B _i which is satisfied in some region of B, i.e. B min < B < B max . The 
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lower limit B min in this region determined by inequality > Mb — Mb-i when the 
system becomes unstable with respect to decay to the nucleons. The upper limit B max 
is determined by the region of applicability of our approach when the baryon density 
in the bulk becomes close to the nuclear density, and therefore, our calculation scheme 
(based on description in terms of quarks ) becomes unjustified at this point. Different 
approaches (based on consideration of hadronic rather than quark degrees of freedom) 
have to be used in this regime. It could happen that some metastable (or even stable) 
states may exist in this low-density regime. However, the corresponding analysis is 
beyond the scope of the present work and it shall not be considered here. 

Therefore, if sufficiently large number of quarks (determined mainly by the axion 
properties) is trapped inside the axion bubble during its shrinking, it may result in 
formation of an absolutely stable QCD-ball with the ground state being a diquark 
condensate. Such QCD-balls, therefore, may serve as the cold dark matter candidate 
which amounts about 30% of the total matter/energy of the Universe, VL DM ~ 0.3 [1J - 

Strictly speaking, the QCD-balls being the baryonic configurations, would behave 
like nonbaryonic dark matter. In particular, QCD-balls, in spite of their QCD origin, 
would not contribute to fish 2 ~ 0.02 in nucleosynthesis calculations because the QCD- 
balls would complete the formation by the time when temperature reaches the relevant 
for nucleosynthesis region T ~ lMeV. Once QCD-balls are formed, their baryon charge 
is accumulated in form of the diquark condensate, rather than in form of free baryons, 
and in such a form the baryon charge is not available for nucleosynthesis. Therefore, the 
observed relation Qb ~ &>dm within an order of magnitude finds its natural explanation 
in this scenario: both contributions to Q originated from the same physics at the same 
instant during the QCD phase transition. As is known, this fact is extremely difficult 
to explain in models that invoke a dark matter candidate not related to baryons. 

Before we continue the description of our proposal we would like to make few 
comments on what have happened on the theoretical side during the last few years, 
which are crucial elements in our present discussions, and which were not available to 
researchers earlier. 

First of all, there existence of the axion domain walls, related to the symmetry 
under discrete rotations of the so-called 9 angle 9 — > 9 + 2im has been known for a long 
time since ^0]- However, the structure of the domain wall considered in had only 
one typical scale, m^ 1 ^> lfermi. Therefore, the quarks, even if they were trapped 
inside the bubble at the very first moment, could easily penetrate through such domain 
wall configuration during the bubble evolution. In this case the axion domain wall ( 
without support of the fermi pressure from the bulk) would completely collapse. What 
was realized only quite recently, is the fact that the axion domain walls have actually 
sandwich substructure on the QCD scale Aq^ d ~ lfermi. Therefore, the fermions 
which are trapped inside the bubble at the very first instant, can not easily penetrate 
through the domain wall due to this QCD scale substructure, and will likely stay in 
the bulk, inside the bubble. In this case, the collapse of the axion domain wall stops 
due to the fermi pressure in the bulk. The arguments ( regarding there existence of 
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the QCD scale substructure inside the axion domain walls) are based on analysis °f 
QCD in the large N c limit with inclusion of the rf fields and independent analysis 
of supersymmetric models where a similar 9 vacuum structure occurs. 

The second important element of our proposal not available earlier, is related 
to the recent advances [3], [I] in understanding of CS phase. The fact that the color 
superconducting phase may exist at high baryon density was discussed a while ago|2j, 
however it was not a widely accepted phenomenon until recent papers [3] where a 
relatively large superconducting gap A ~ lOOMeV with a large critical temperature 
T c ~ 0.6A were advocated. 

To conclude the Introduction we should remark here that the idea that some quark 
matter, such as strange quark "nuggets" may play a role of the dark matter, was 
suggested long ago [13], see also original papers [Hj and relatively recent review [T5j on 
the subject. The idea that soliton-like configurations may serve as a dark matter, is also 
not a new idea[in|. Most noticeable example is being Q-balls[T7]. The idea that the dark 
matter may be just solitons containing large baryon (or even antibaryon) charge is, again, 
an old idea^H], see also ^H]. The new element of this proposal is the observation that 
one can accommodate all the nice properties (discussed previously [EJ- [TH]) without 
invoking any new fields and particles (apart from the axion). Rather, our QCD-balls 
formed from the ordinary quarks which however are not in the "normal" hadronic 
phase, but rather in color superconducting phase when squeezed quarks organize a 
single coherent state described by the diquark Bose-condensate, similar to the Cooper 
pair condensate in BCS theory in conventional superconductors. 

In many respects ( in terms of phenomenology) the QCD balls are similar to 
strangelets|13j-[T5] with few important differences, see below: 

• 1. In our proposal the first order QCD phase transition is not required for the formation 
of the QCD-balls. Axion domain walls of a large size (in comparison with a typical QCD 
scale) are able to form the large bubbles. These bubbles, filled by u, d, s quarks, play 
the same role as the bubbles formed during the first order phase transition as discussed 
in[IBJ. 

•2. The Stability of strange quark matter at zero external pressure, as described in|13j- 
[T5] . is highly model dependent result. In particular, the stability of strangelets is very 
sensitive to the magnitude of the bag constant within MIT bag model calculations. The 
idea which is advocated in the present work has a new element, the external pressure 
due to the axion domain walls. With this new element the stability of the system is very 
likely to occur in very wide region of the parametric space even in the models which 
would not support strangelets in the absence of the external pressure. 
•3. The bulk of the QCD ball is in the superconducting phase. This property obviously 
influences the phenomenology of how the QCD balls interact with a normal matter. 
In particular, if the energy of a hadron which hits the QCD ball is smaller than the 

% Uniqueness of the rj' field in this problem is related to the special structure of interaction of the axion 
held 6{x) and the singlet rf{x) held in low energy description of QCD when only a special combination 
[9(x) ~ rf{x)\ is allowed to enter the low energy QCD Lagrangian. 
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superconducting gap A, the hadron can not penetrate into the bulk and excite the 
internal degrees of freedom of the system, but rather it will be reflected (the so-called, 
Andreev reflection). Similar property is also true for the strangelets |1 Hj . The elastic 
cross -section of hadrons on the QCD balls is large, of the order of the geometrical size 
of QCD balls; the inelastic cross -section (when internal degrees of freedom are excited) 
is almost identically zero for small energies as mentioned above. Electromagnetic 
interactions of photons with QCD balls contain the standard fine structure constant a 
with an addition suppression due to the neutrality of the CFL dense quark matter. Such 
features for the interactions imply that if the QCD ball with small velocity v/c ~ 10 -3 
enters the Earth, it will not decay by exploding. Rather it will go through the Earth and 
exit on the opposite side of the Earth leaving behind the shock waves. It is tempting to 
interpret the recent seismic event with epilinear source [20] as the process which involves 
the dark matter particle, similar to the QCD ball. 

•4. There is a maximum size of the QCD-ball above which such an object can not be 
formed and can not be absolutely stable. This is due to the fact that for very large system 
the axion domain wall pressure becomes a negligible factor which can not stabilize the 
system. 

•5. The property on a maximum size mentioned above has a profound phenomenological 
consequence. Indeed, if one assumes that stable state as described in[l3] exists, the 
strangelets can collide with an ordinary neutron star which results in formation of a 
quark star. In such a case all neutron starts would be transformed into quark stars. 
In our proposal, when the maximal size of the QCD ball is determined by the external 
axion domain wall pressure, this transition (from neutron stars to quark stars) does not 
happen routine effect. 

•6. If the size of the QCD ball slightly exceeds the maximum critical size, it becomes 
metastable, rather than stable configuration. Such QCD balls could also be interesting 
particles for the dark matter phenomenology, see footnote on page 14 and discussions 
in the Conclusion. 

2. QCD-balls 

Crucial for our scenario is the existence of a squeezer, axion domain wall which will 
be formed during the QCD phase transition. As is known, there are many types of 
the axion domain walls, depending on a model. We assume that the standard problem 
of the domain wall dominance is resolved in some way as discussed previously in the 
literature, see e.g. [21], and we do not address this problem in the present paper§. We 
also assume that the probability of formation of a closed bubble made from the axion 

§ It is widely accepted that the domain walls in the so-called, N=l axion model will be eaten up by 
the axion strings at a very high rate. That is true for the axion walls bounded by strings. However, 
if a domain wall is formed as a closed surface, the probability for such a wall to decay is extremely 
small. Therefore, such domain walls in N = 1 model can play the same role in our scenario as stable 
domain walls in N ^ 1 models. Besides that, N = 1 model has a nice property that the domain wall 
dominance problem is automatically resolved. 
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domain wall is non-zero||. We also assume that quarks which are trapped in the bulk, 
can not easily escape the interior when the bubble is shrinking. In different words, the 
axion domain wall is not transparent due to the QCD sandwich structure of the wall as 
discussed in |H],[I21- The collapse is halted due to the Fermi pressure. Therefore, we 
assume that a large number of quarks remains in the bulk, inside the bubble when the 
system reaches the equilibrium. 

2.1. Equilibrium 

The equilibrium is reached when the Fermi pressure equals the surface tension and 
pressure due to the bag constant E B . To put this condition on the quantitative level, 
we represent the total energy E of a QCD-ball with the fixed baryon charge B, in the 
following way, 

E = 4naR 2 + ^R 3 + ^E B R 3 (1) 

07T 3 

where we assume the quarks to be massless. We assume that the relativistic fermi gas 
is non-interacting in the first approximation, see corrections due to the interactions 
below. In this formula fi is the Fermi momentum of the system to be expressed in 
terms of the fixed baryon charge B trapped in the bulk; R is the size of the sysytem; 
g is the degeneracy factor, g ~ 2N c Nf = 18 for massless degrees of freedom; E B is 
bag constant which describes the difference in vacuum energy between the interior and 
exterior. The bag constant is a phenomenological way to simulate the confinement. 
Finally, a ~ fa^nfn is the axion domain wall tension with f a ~ (10 10 — 10 12 )CreV being 
constrained by the axion search experiment. 

In what follows, it is convenient to introduce dimensionless scaling variable x, as 
follows, x\f~B = R\[E~b such that energy per quark etot = E/B can be expressed in the 
following simple way in terms of dimensionless parameters x and o~o, 

, . E ^1/4 ( A n 3 ,/97T 47T ,\ ,„ s 

£ „ W ^ = 4 /4 (W + -^- + T ^J (2) 

E 1 ^ 

X = R WJ^ a ° = B^Ef- 

The minimization of this expression de to t{x) /dx\ x=xo = determines the stability radius 
Xq which fixes the energy of the system at the equilibrium, e ta t(xo). In particular, if one 
neglects o"o in eq. (J2J) originated from the axion domain wall tension, one reproduces 
the well known results, x ~ 0.48, e(x ) ~ 1.9E B 4 . Such a relation means that if E B 
is relatively small such that the energy per quark is less than mAr/3, the configuration 
becomes an absolutely stable state of matterfTSj-fTK]. 

| We do not attempt to develop a quantitative theory of the formation of the QCD-balls in this work; 
It is sufficient for our following discussions that this probability is finite. 
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In eqs. JJJ [2]) we have neglected many important contributions which can drastically 
change the results. We shall review the role of these contributions below. The main 
goal of this subsection is the incorporation of these contributions into eqs.((TJ Ej). First 
of all, in eq. (0) we neglected the quark-quark interaction on the Fermi surface, 
which brings the system into superconducting phase for relatively large baryon density 
[3]. The corresponding contribution AE int to the total energy ((TJ) is negative and at 
asymptoticaly large fi is equal to| 



AE int = f- ■ (-R 3 ) (3) 

The negative sign of AEi nt is quite obvious: the formation of the diquark condensate 
due to the quark-quark interaction lowers the energy of the system. For appropriate 
treatment of this term one should express fi as a function of B, R according to the 
relation ((TJ) and substitute this into eq. ©• In principle, one should also take into 
account that the superconducting gap A(/i) also strongly varies with \i (and therefore, 
with R) in the relevant region of fi. However, in what follows we shall ignore this 
dependence for numerical estimates and shall treat A ~ lOOMeV as constant. Our last 
remark regarding eq. (JBJ). This formula was derived for very large fi. Nevertheless for 
illustrative purposes we shall use the expression for AE int for small /i as well. We shall 
see that in the relevant region of densities the contribution AEi nt does not exceed 15%. 
This somewhat justifies the use of expression (J3J) for our numerical estimates which 
follow. With all these reservations in mind, we account the additional contribution to 
energy per quark, describing the quark-quark interaction on the Fermi surface by adding 
Ae 1 ^ into eq. (J2J) in the following way 

Ae fot = -Ell 1 | \iz ' ~Tr= ••'')■ (» 




where we expressed everything in terms of dimensionless parameter -7^= and 
dimensionless variable x. 

Now we want to consider the modification of eq. (0) which is related to the actual 
variation of the bag "constant" Eb with fi. To explain the physical meaning of this 
effect, we remind the reader that the bag "constant" Eb describes the differences of 
vacuum energies in the interior and exterior regions. It is a phenomenological way to 
simulate the confinement. The bag "constant" contribution goes with the positive sign 
to E, see eq.((TJ). The physical reason for this sign is obvious: the vacuum energy outside 
the bubble is lower than inside, thus the positive contribution to E, in contrast with 
the interaction term, — — £ discussed above. 

Our main point is as follows: the contribution related to Eb can be expressed 
formally in terms of the difference between the vacuum condensates calculated at zero 
(exterior) and non-zero (interior) baryon densities. The most important contribution to 
Eb is due to the gluon condensate, such that Eb{(j) ~ (1§|^/L)m=o — (tll^V^A^o with 
b = yiV c — |iVj where we used the well-known expression for the conformal anomaly 
in QCD in the chiral limit. We do not know Eb{ij) as a function of fi for the relevant 
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region of the baryon density. However we do know the behavior of this quantity for 
relatively small densities corresponding to the nuclear matter densities [~ 



(fgLWo (O.Q5GeV)p N p N 

QG&r* " (v^),=o " {2UMeVf [ } 

where p^ is baryon density, and the magnitude for the gluon condensate is known to 
be, {^G 2 )^ = q ~ 1.2 ■ 10~ 2 GeV A . As expected the gluon condensate (and therefore, 
the absolute value of the vacuum energy) decreases when the baryon density increase. 
Similar formulae are known for the chiral quark condensate where for the small densities 
one can derive the following relation ^_ q \^° = 1 — aN 2 p S with sigma term measured to be 
<tat ~ A5MeV see j2Hj for the details. One should emphasize here that the formula 
describing the variation of the gluon vacuum condensate at small baryon densities pjy, 
is a direct consequence of the QCD low energy theorems. It is a firm result of QCD, 
not based on any model dependent considerations, and should be accepted as it is. 

More specific information on the bag "constant" Eb contribution as function of 
p in the entire region of of p can be calculated in some non-physical models such as 
QCD with two colors, N c = 2 24 . Such a knowledge can not be literally used for our 
numerical estimates which follow, however it can be used for modeling the functional 
dependence of the vacuum energy. 

Therefore, we want to model two properties discussed above in order to incorporate 
them into the corresponding eq. First, the bag constant contribution must vanish 
when the baryon density in the bulk vanishes. This corresponds to the case when 
vacuum energy inside and outside of the bubble is the same, and therefore, it should be 
no additional vacuum energy contribution to the equation for the equilibrium. Secondly, 
the bag constant contribution should vary with density as we discussed above. 

Our first parametrization is motivated by analysis t 24j of the vacuum condensates 
in QCD-like theories at finite baryon density as a function of p. If we assume a similar 
behavior in real QCD than we should replace the bag constant Eb by the expression 

2 

E B — > E B (l — %) for p > p c and Eb — > for p < p c , where p c would correspond to 
a magnitude of the critical chemical potential at which the baryon density vanishes. In 
QCD, one expects that this is to happen at p c ~ 330MeV. 

As before, one should express the corresponding contribution to e to t m terms of 
fixed baryon charge B and radius R, such that the bag" constant" contribution actually 
becomes a complicated function of B,R. In terms of dimensional parameter x the 
corresponding contribution to (J2J) is accounted for by the following replacement, 

4 /4 -* 3 =► ET-x* .(l- (V 3 ■ J^=A (6) 



o o \ vi y Eb i 

Let us emphasize: we are not attempting to solve a difficult problem of evaluation of 
nonperturbative vacuum energy as a function of p in QCD. Rather, we want to make 
some simple estimates to account for this effect in order to analyze the stability of QCD 
balls later in the text. 

We want to be confident that the results on stability of QCD balls (to be discussed 
later) are not sensitive to the specific parameterization © motivated by the study 
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of QCD with two colors. Therefore, we would like to have a different, independent 
parameterization of the same effect to be used in our stability analysis. We make use 
of eq.(JHJ) which is valid for small densities pn- This formula gives us an idea about 
typical variation of vacuum condensates when the baryon density changes. We assume 
that the vacuum energy difference in QCD (the bag "constant" contribution in eq. (J2J)) 
can be expressed in terms of different vacuum condensates with the typical scale for the 
variation given by eq.©. 

We want to implement the QCD property © into the MIT bag model. If the 
phenomenological numerical magnitude for the bag constant Eb were closed to the 
numerical value for the vacuum energy (ff^Gj \ v ) ~ (340Mey) 4 we could literally 
use eq ©, such that the bag constant contribution can be parameterized as follows, 
Eb(pn) — Eb (264MeV r ') 3 • Unfortunately, these two are very different numerically, 



and we will introduce the corresponding correction factor r = ^(^G^/Eb — 
(3A0MeV) I \l50MeV) ~ 2.25 in our implementation of QCD property (© into the 
MIT bag model, see below. 

Still, formula Eb(pn) ~ Pn can not be used literally for our purposes because we 
need an expression for the bag "constant" contribution which goes to constant Eb at 
large densities, E b {pn) —> E B - A simple model which satisfies this requirement is to 
make the following replacement, 



where we introduced the correction factor r to match the scales. As before, one 
should express the bag "constant" contribution proportional to ((Zj) in terms of a fixed 
baryon charge B and radius R. We shall analyse the corresponding equation (j2J) with 
improvements (|7j) in the next subsection. To anticipate the events, one should mention 
that our two models (JB1 Ej) describing the effect of the bag "constant" variation with 
baryon density lead to the similar results, see below. 

The next approximation we have made in eqs.(P][2J) is related to the assumption of a 
thin-wall approximation for the domain wall. This may not be well justified assumption 
because the typical width of the domain wall and the size of QCD ball could be the same 
order of magnitude, such that thin-wall approximation is failed. However, we neglect 
these complications at this initial stage of study. Nevertheless, we do not expect that 
this effect can drastically change our qualitative results which follow. 

We also neglected in eqs.(HJ 121) ah complications related to the finite magnitude 
of the quark masses, first of all m s , which result in additional K condensation along 
with diquark condensation in CFL phase (2Hj. Finally, the expression for the energy E 
with corrections (|3I6|) . changes the simple relation (© between baryon charge B and 
chemical potential p according to the standard thermodynamical relations, B — — 
where F = E — pB is the free energy. However, we checked that these changes are 
relatively small( do not exceed 5% in the relevant region of p). Therefore, in what 
follows, in order to avoid the technical complications in the qualitative analysis, we 





r 3 p N 



Eb 



(7) 
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use a simple algebraic expression ([TJ which is formally valid only for noninteracting 
quarks, B ~ fi 3 , but numerically remains a good approximation in a large region of 
/i. This allows us to use the dimensional variable x which we introduced before for 
the non-interacting case. Let us repeat again: we do not attempt to solve the problem 
quantitatively with all uncertainties in parameters discussed above; rather, we want to 
give some qualitative arguments demonstrating that stability region might occur in the 
wide region of B with realistic choice of parameters specified below. 

With all these reservations regarding eqs.^EJ) in mind we express the energy of a 
QCD-ball per baryon charge B in units of \/Eb , as follows 

y(x) tot = E^eUx) = ^\l- {-f' 3 ^-x> ] (8) 




+ 1 Ana ° x2 + h 

In this formula, in comparison with eq.(J2J), we took into account the effect describing 
the quark-quark interaction on the Fermi surface given by eq. (J3| and the effect of the 
variation of the vacuum energy with baryon density, given by eq. (jHJ). 

The equilibrium condition de to t(x = x )/dx = determines the radius x of the 
QCD ball with baryon charge B. We shall analyze this condition in the next subsection; 
now we want to constraint Xq < x to be considered. The constraint follows from the 
condition that the baryon density should be relatively large. In this case our treatment of 
the problem by using the quark degrees of freedom, eq.(|HJ), rather than hadronic degrees 
of freedom, is justified. The baryon number density p^ for the QCD ball configuration 
is given by^f, 

B E 3/4 

which gives upper limit x above which our approach is not justified. Numerically, with 
our choice of parameters, see below, x ~ 0.6, and therefore, any solution xq of the 
equilibrium condition de tot (x = Xq)/8x = must satisfy to the constraint xq < x ~ 0.6. 



2.2. Stability of QCD balls 

As expected, the equation describing the equilibrium de tot (x = x )/dx = has a 
nontrivial solution (minimum) in a large region of parametrical space deterimed by 
parameters Eb, ct, A, /i c , B. It is not our goal to have a complete analysis of this allowed 
region of solutions. Rather, we shall make a specific choice for all parameters except 
for the baryon number B and analize the stability condition as a function of B. We 
shall also comment on results with a = corresponding to pure QCD configuration 
without any involvement of the axion field (case considered previously in MIT bag 
model, [13]- [IS])- The first step is to calculate the point x = xq which is determined 

Our normalization for the baryon charge corresponds to B = 1 for the quark, thus factor B/3 in eq. 
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by equation de tot (x = x )/dx = 0. The next step is to analyze the stability of the 
obtained configuration as a function of external parameters. Condition when the QCD- 
ball becomes an absolutely stable object can be derived from the following arguments. 
Total energy per quark e tot (x ) in eqs. (0 E|) is a combination of two factors: the 
first one, €qcd(xo), is due to the strong interactions; the second factor, e ax i on (xo) is 
mainly due to the axion domain wall tension" 1 ", i.e. e to t(x ) = €qcd(%o) + £axion(%o) , with 
£ a xion( x o) = (^fo^o) an d £qcd(xq) is determined by rest of terms in eq. (JHJ). The 
absolute stability of the system implies that a nucleon can not leave a system because 
the energy of the configuration with baryon charge B is smaller than the energy of 
configuration of charge B — 3 plus energy of a nucleon with baryon charge B = 3 and 
energy of the axion emission. Such a situation is analogous to the three dimensional 
quantum mechanical problem with an effective potential being a step-function and the 
energy of the bound state is lower than the potential energy at the large distances. In 
this particle obviously can not leave the system. 

We should emphasize here that the quarks can not leave the system due to the 
energetic conditions which take place after the QCD ball is formed. In different words, 
the stability occurs due to the differences in properties inside/outside of the QCD ball, 
and not due to the features of the original axion domain wall. The axion domain wall 
already had played its role during the formation period when a large number of quarks 
could not escape the system and were trapped in the bulk during the collapse of the wall. 
A similar situation when a configuration may become a stable one due to a difference 
in conditions (inside/outside the bulk) was discussed long ago [IE] as an example of a 
non-topological soliton in quantum field theory. We further comment on the similarities 
with non-topological solitons later in the text. 

It is quite obvious that the axion domain wall with a typical correlation length 
~ m" 1 ^> Aq 1 cd can not produce nucleons by itself when it shrinks due to the nucleon 
emission. Instead, typically, the axion domain wall reduces its size by emitting the axions 
while the nucleon leaves the system. In this case the term e axion (x ) = E l ^ A (Aira xl) 
is responsible for the emission of axions rather than production of nucleons. As a 
result of this, this term should be ignored for the analysis of the stability. However, 
with exceedingly small probability the emitted axion, in principle, can be absorbed 
by the nucleon which leaves the system. In this case the energy, in principle, can 
be transformed from the axion domain wall to the produced nucleon and the term 
£axion{xo) = E^ 4 (4naoxl) should be accounted in the energy budget for analysis of the 
decay. We estimate in appendix that the probability for the corresponding absorption 
of the axion by the leaving nucleon is negligible. Therefore in what follows we neglect 
this process. 

The relevant term which describes the emission of nucleons is the one related to 
the QCD physics i.e. €qcd(xo)- Therefore, the condition when configuration becomes 
a sufficiently stable (with the life time exceeding the life time of the Universe, see 

+ the QCD contribution to a due to the rj' and pions is suppressed by a factor f%/ f% <C 1. 
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Appendix for details) is determined from the following inequality 

m N de tot {x) _ 

tQCD{XQ) < — , 7^ \x=xo =0, x < x, (10) 

where the last condition follows from 

To analyse eq. (|10|) we shall accept the following magnitudes for the dimensional 
parameters: 

A^lOOMel/; a ~ 1.8 • 10 s GeV 3 ; (11) 
// c ~ 330Me\/; E B ~ (150Me\/) 4 . 

Having these external parameters fixed, we left with the only one unknown number, the 
baryon charge B, which eneters ctq in our dimensionless parametrization (|2|8)l . We shall 
treat do as a free parameter and our goal is to find the region of a when conditions 

are satisfied. As we discussed above, we shall use two different models to account 
the effect of the variation of the bag constant contribution with density, see eqs. (®[7J). 

Having defined our stability condition ([TO]), external parameters ([TTjl and two simple 
models accounting the effect of the variation of the bag constant, eqs. (®IZJ), we reduce 
our problem to analysis of dimensionless functions, Uqc D {x) and yQc D (x) defined as 
follows, see eqs. (jEHHIEl), 

ytot(x) = yQc D (x) + y aX ion(x); y a xion(x) = 4:iro- x 2 (12) 
Vqcd = + ^YT^ ~ °' 48x ' (13) 

Vqcd = ^— + 4.2a; 3 (l - 5.68x 2 ) - 0.48x, (14) 

oc 

where three consequent terms describe: the fermi pressure, the bag constant contribution 
accounting the variation of the vacuum energy with the baryon density ()7I6|) . and, finally, 
the quark-quark interaction on the fermi surface (@J) correspondingly. Stability condition 
(|10|) in dimensionless variables becomes 

(i,2)/ \ ^ ™n 01 dy t Qt(x) 
yqcb^o) < ^= ^ 2.1, — \ X=X0 = 0. (15) 

Before we discuss some specific numerical results which follow from analysis of eqs. (fTBI 
JTHj) . we would like to list some general model-independent properties of the solutions. 
We believe that the properties listed below are quite common features of the QCD balls, 
which likely to remain untouched even in a more general treatment of the problem when 
many additional effects are included (some of these effects were mentioned above). 

a). As we already mentioned, in the absence of the axion field, a = 0, the problem 
was extensively discussed earlier using MIT bag model. [T3* ] -|15j. Our original remark 
here is: when a variation of the vacuum energy with density is taken into account, a 
stable solution disappears provided that a typical QCD scale for the vacuum variation 
()6I7|) is used. The physical reason for that behavior is quite obvious: a density - 
dependent vacuum energy is not a sufficiently strong squeezer to equilibrate the fermi 
pressure. A typical scale for the variation should be reduced (in comparison with what 
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Legend 



mu=5.68 

mu=0.8 

mu=0.3 



Figure 1. In this figure we plot the value of the total energy of the QCD ball per 
baryon charge I|12I14|I with a = 0. Three curves corespond to the different parameters 
(■|) 2 / 3 -^£= = 5.68 ,0.8 ,0.3 modeling the variation of the bag constant with the 
baryon density. The main observation: a minimum corresponding to the equlibrium 
does not exist for the physical parameters when (^) 2 / 3 -^= = 5.68 from eq. I|14|l . 
Equlibrium appears when a typical scale for the variation is reduced by an order of 
magnitude. 



we assumed in eqs. (l6!7j) ) by an order of magnitude, in order for the solution to reappear. 
Specifically, we checked that the equilibrium is possible for a = if coefficient 5.68 in 
(I14J) describing the vacuum energy variation is replaced by 0.5. We demonstrate this 
effect in Fig.l where we display the total energy of QCD ball per baryon charge (fT2|) with 
(j = for three different values of parameter (| ) 2 ^ 3 ^g^ describing the effect of variation 
of vacuum energy with the baryon density. For a typical choice of physical values (jllj) 
the relevant parameter is (^) 2//3 "^^ = 5.68. In this case the minimum does not exist 
which implies that the stability can not be achieved as announced above. The minimum 
starts to reappear only when a typical scale for the variation of density is considerably 
reduced, see Fig.l with the curve corresponding (~) 2//3 "^7g^ = 0.3. Such a small value for 
the critical value fj% does not look appealing from the physics point of view. Therefore, 
we incline to accept that there is no solution for such a configuration (strange quark 
nuggets. |13j-|15j) in QCD if no external pressure (such as gravity or axion domain wall) 
is applied. It is certainly not a very new result: special study on stranglets reveals[26 
a strong model dependence of the stability of strange quark matter. In particular, the 
Nambu Jona- Lasinio model does not support any kind of strangelets[2Ij. 

b). In general, one expects there existence of a minimal and maximal sizes for 
the QCD balls in the region of stability. The minimal charge B m i n corresponds to the 

1 /3 /fTTTl. 

maximum a™ ax ~ B min . At this point the the stability requirement (jlOJ) is marginally 
satisfied. When B < B min , ctq becomes too large such that nucleons can leave the 
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Figure 2. In this figure we plot the value of the total energy of QCD ball per baryon 
charge (|12I13|I for Anao — 1,5,10. The main observation: a minimum describing the 
equilibrium at Attctq = 1 corresponds to the maximum possible baryon charge. Solution 
goes away for smaller a®. The equilibrium at Airao = 10 corresponds to the minimum 
possible baryon charge B m i n when solution at the equilibrium still satisfies the stability 
requirement (|10|l . At larger <tq > <r™ ax the quark energy per baryon charge becomes 
large enough such that nucleons can leave the system. 



system. On the other hand, the maximum possible charge, B maxj corresponds to the 
minimum value of <r™ n ~ B^]/ x 3 . For larger B, the baryon density Q becomes too 
low to justify our approach (based on the quark degrees of freedom). At lower baryon 
densities some metastable states may form; they could decay to some heavy elements 
which might be of interests for astrophysics. However the corresponding study would 
require an analysis of the system in terms of nuclear degrees of freedom, which is beyond 
the scope of the present work. When <Jo becomes even smaller, the problem is essentially 
equivalent to a = studied earlier where stable solutions are not expected to occur. 

Numerically, we analyzed two models (fTT^ IT4*j) which lead to the similar results. In 
particular, for model ()13|) the maximum possible tension, 47ra™ ax ~ 10 corresponds 
to the minimum baryon charge B m i n . For such do the equilibrium is reached at 
x — 0.32 wnen the energy per quark i/qcd( x o) — 2.1 hits the upper energy bound 
of the stability region (fTK|) . When 47ra™ ax > 10, the energy per quark becomes too high 
such that nucleon can escape and the system would decay. In physical units this solution 
corresponds to B min ~ 10 32 and stabilization radius Rq = xq\[B j \[E~b — l0 11 GeV~ 1 . 



Energy per quark for this configuration eq CD = 
is smaller than constituent quark mass, as it should be*. 

For the same model, the minimum possible tension when our approach is justified, 
47rcr™ m ~ 2 corresponds to the maximum possible baryon charge B max . According to the 

* We remind that we discuss the QCD part of energy only; the total energy of the configuration which 
includes the axion part is larger. 
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scaling B ~ cr 3 ,, the maximum baryon charge B max = (^itr) 3 B min ~ 10 34 is two orders 
of magnitude larger than B m i n . In this case the equilibrium is reached at xq — 0.52 
when the baryon density (JUJ) is already relatively low, and close to the boundary when 
the quark based lore can not be trusted. 

Our second model $T^i gives quantitatively similar results, and it is not worthwhile 
to discuss numerical details here. The most important features of the solution for this 
model remain the same: there is a region between B max and B min when solutions are 
stable; at o = solution does not exist at all provided that a typical QCD scale for the 
vacuum variation (j6!7J) is used. 

However, one should take all these numerical estimates very cautiously because of 
a number approximations we have made in eqs. Nevertheless, in what follows, 

mainly for the illustrative purposes, we shall stick with these numerical estimates. 

The quark number density n in the region B m i n < B < B max when our approach 
is justified is estimated as 



where we used the expression (0) for the baryon density. As we already mentioned the 
expression (^) is formally valid only for noninteracting quarks, but numerically remains 
a good approximation in a large region of /z. In eq. (J 16)) 3n ~ 3(108MeV A ) 3 , is the 
nuclear saturation density normalized with our convention ( B = 1 for quarks), thus 
factor 3 in front of the numerica value 0.16(/m)~ 3 ~ (lOSMeV") 3 . It is quite remarkable 
that the numerical value for n is in the region where color superconductivity phase is 
likely to realize, and therefore, our treatment of the squeezed fermi system as the quark 
dense matter (rather than ordinary nuclear matter) is justified a posteriori. 

Few remarks are in order regarding eq. lfTo^) . First of all, the estimates presented 
above demonstrate that we are in the region of the phase diagram where CFL phase 
is likeley to realize. Therefore, our original assumption is justified. Secondly, for large 
B > B max our treatment of the system is not valid anymore, and a different type of 
QCD balls with an ordinay nuclear matter (instead of diquark condensate) in the bulk 
may be formed and could be even stable in some regions of parametrical space. Though 
this region of large B > B max could be an interesting region from the phenomeno logical 
point of view, it shall not be discussed herejj. However, even in this case when the QCD 
balls made of nuclear matter, rather than quark dense matter, we still expect that there 
should exist a maximum size above which the stability is not possible. This follows from 
our analysis that stability can not be achieved without the external pressure P a due to 
the axion domain P a ~ 2a /R which vanishes at very large R. 

Another factor which also constraints the size of the balls is related to the 
suppression of large size closed axion domain walls during the formation stage. It is 
clear that the formation of the large size closed domain walls is suppressed according 

jj The corresponding proper treatment would require the knowledge of the dynamics of the interacting 
nuclear matter, which is not the subject of the present work. In principle such nuclear matter could be 
also stable. 




(16) 
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to the Kibble-Zurek mechanism [25], |2Hj; however an explicit estimation for this effect is 
still missing. 

As we mentioned in the Introduction, we do not address the problem of formation 
of QCD balls in this letter, it will be a subject of a different work. However we would like 
to mention some relevant elements of a possible scenario of how QCD-balls, in principle, 
can be formed after the QCD phase transition, at a temperature of order 150 MeV 
which is much higher than the critical temperature for quark pairing estimated to be 
~ 0.6A. The main point is this: the axion domain wall with the QCD-scale substructure 
as discussed in[ll is very selective with respect to the momentum of the particles; it 
is almost transparent for light it mesons with large momentum k > m„ such that the 
transmission coefficient is close to one. Therefore, the highly energetic pions can easily 
penetrate through the domain wall and leave the system. 

At the same time, the transmission coefficient is close to zero for slow-moving 
particles such as baryons with k < m n . Eventually, this "selective" feature of the domain 
wall may cool down the system considerably. Due to the domain wall pressure it may 
reach the critical density when it undergoes a phase transition to a color superconducting 
phase with the ground state being the quark condensate. At this point we assume that 
the baryon number trapped in the bulk is sufficiently large. If B 3> B max , the quarks 
will leave the system by forming nucleons until the upper limit B max is achieved. At 
this point the energy per unit baryon charge €qcd = v' ' EbIIqcd{.xo) < becomes 
sufficiently small such that quarks can not leave the system. Some specific calculations 
are required before any statements regarding a possibility to form the QCD balls can 
be made. We do not see any fundamental obstacles which would prevent the formation 
of such objects. Terefore, at this moment we simply assume that this probability does 
not vanish. 

2.3. QCD- balls versus Q-balls 

In this subsection we would like to mention a striking resemblance of the QCD-balls 
(which is the subject of this letter) and non-topological solitons[H)], as well as Q-balls [T7] 
which is a special case of a nontopological soliton configuration associated with some 
conserved global Q charge. Both cases, QCD balls and Q-balls demonstrate a similar 
behavior for a soliton mass as function of Q. Namely, QCD balls as well as Q-balls may 
become very stable configurations for relatively large Q charge. Therefore, an effective 
scalar field theory with some specific constraint on the effective potential (when Q ball 
solution exists) is realized for QCD in high density regime by formation of the diquark 
scalar condensate which plays the role of the effective scalar field. The big difference, 
of course, that underlying theory for QCD-balls is well known, it is QCD with no free 
parameters. This is in contrast with the theory of Q-balls when the underlying theory is 
not known. Formal similarity becomes even more striking if one takes into account that 
the ground state of the CFL phase in QCD is determined by the diquark condensate 
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with the following time dependence ~ e l2/rf , 

<*£*S>' ~ K^)' ~ • ^ c e o6c , (17) 

with being the original QCD quark fields, and \i being the chemical potential of 
the system, see formula (40) from ref. |3U| . As is known, such time-dependent phase 
is the starting point in construction of the Q balls [T7]. In the expression (JT7J) we 
explicitly show the structure for the diquark condensate corresponding to CFL (color- 
flavor locking) phasejl] with (a, /3, etc.) to be flavor, (a, b, etc.) color and (i, j, etc.) 
spinor indices correspondingly. Of course, there are many differences in phenomenology 
between Q balls|17j and QCD-balls. For example, in CFL phase the baryon symmetry is 
spontaneously broken, and corresponding Goldstone massless boson carries the baryon 
charge. However, the evaporation of this massless particle into hadronic phase from 
the surface of the QCD-ball is not possible, because hadronic phase does not support 
such excitation. This is in contrast with phenomenology of Q-balls, where the theory 
is formulated in terms of one and the same scalar field, such that evaporation of 
particles from the surface of the Q-ball is possible if some conditions are met. In spite of 
many differences, the analogy with Q-balls is quite useful and can be used for analysis 
of different experimental bounds on QCD-balls, which is the subject of the next section. 

3. Experimental bounds on masses and fluxes of QCD-balls 

In this section we adopt the results of paper [HI] to constraint the free parameter (charge 
B) of the QCD-balls. In the paper [31] the authors re-analyzed the results of various 
experiments, originally not designed for the Q-ball searches, but nevertheless these 
experimental results were successfully used in jH5 to bound different properties of the 
Q-balls. We actually repeat this analysis for a specific type of the QCD-balls when 
original quarks are in the CFL (color-flavor locking) phase[4j. 

As we mentioned earlier, at sufficiently large baryon density, the color 
superconductivity phenomenon takes place. However, there are many different phases 
(as a function of parameters like m s , number of light flavors, etc.) associated with 
color superconductivity. In particular, for 3 degenerate flavors of light quarks, the CFL 
phase with nonzero value for the diquark condensate (fTTj) is realized. Due to the fact 
that equal numbers of u, d, s quarks condensed in the system, the electric charge of the 
ground state is zero, i.e. no electrons required to neutralize the system. This is quite 
important feature for the phenomenology of the QCD-balls we about to discuss. Nature 
is less symmetric, and other CS phases could be realized. In particular, for relatively 
large m s , along with diquark condensate, the K condensate may also be formed[25j. In 
the limit of very large m s , QCD becomes effectively a theory with two light quarks. In 
this case, the Cooper pairs are ud — du flavor singlets. This phase, the so-called 2SC 
(2 flavor super-conductor ) phase is a phase with non-zero electric charge. Electrons 
neutralize the system, however, all properties, such as interaction cross sections, the rate 
of energy loss of QCD balls in matter, are very different for QCD-balls with quarks in 
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CFL or 2SC phase. In what follows, to avoid many complications, we limit ourself with 
analysis of QCD balls where quarks are in the most symmetric CFL phase, in which 
case the QCD-ball has zero electric charge. 

We assume, in analogy with[3T], that a typical cross section of a neutral QCD- 



ball with matter is determined by their geometrical size, ttRq. In this case, the only 
information we need to constraint the QCD-ball parameters, is its size and mass. We 
also assume that the Q CD-balls is the main contributor toward the dark matter in the 
Galaxy. Their flux F then should satisfy 

F < F DM ~ ^ ~ 7.2 ■ 10 5 ^cm- 2 S ec-V-\ (18) 
4:71 M b Mb 

where pdm is the energy density of the dark matter in the Galaxy, p DM ~ 0,3 l v , and 
v ~ 3 • 10~ 3 c is the Virial velocity of the QCD-ball. We identify Mb in the expression 
()18|) with the total energy E of the QCD ball at rest with given baryon charge B. 
The Gyrlyanda experiments at Lake Baikal reported that the flux of neutral soliton-like 
objects has the bound jH2] 

F < 3.9 • 10~ 16 cm" 2 sec' 1 sr~\ (19) 

which translates to the following lower limit of the neutral QCD-ball mass Mb and 
baryon charge B, 

M exp > 2 . 10 21 GeV - ( 20 ) 
S -P^(^£)9/8[^( 87rc 2)i]-9/8 > L6 .10 2 °. 

a v ' A 2 

Similar constraints follow from the analysis of the Baksan experiment |33j and analysis 
[3*T] of the Kamiokande Cherenkov detector [HI], and we do not explicitly quote these 
results. These experimental bounds are well below the critical line of the stability of 
the QCD-balls. 



4. Discussions and Future directions 



Complete theory of formation of the QCD-balls is still lacking. Only such a theory 
would predict whether QCD-balls can be formed in sufficient number to become the 
dark matter. Such a theory of formation of the QCD balls would answer on questions 
like this: 1. What is the probability to form a closed axion domain wall with size 
£ during the QCD phase transition? 2. How many quarks are trapped inside the 
domain wall at the first instant? 3. How many quarks will leave the system and how 
many of them will stay inside the system while the bubble is shrinking? 4. What is 
the dependence of relevant parameters such as: size baryon number density n(t) 
and internal temperature T(t) as function of time? 5. Do these parameters fall into 
appropriate region of the QCD phase diagram where the color superconductivity takes 
place? 6. What is the final density distribution of the QCD-balls as a function of 
their size R after the formation period is complete? 7. Will the QCD balls survive the 
evaporation and boiling even if they formed? Clearly, we do not have answers on these, 
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and many other important questions at the moment. All these interesting, but difficult 
questions are obviously beyond the scope of the present work, and shall not be discussed 
here. However, we want to make a short comment on issue 7 which was an important 
element in many previous studies. 

The question on evaporation of quark nuggets was discussed earlier, see original 
papers [US]- [HE] and recent review [T5| . The first study of this question is due to Alcock 
and Farhi[3S] who argued that only very large nuggets with B > 10 52 could survive 
the evaporation. This result would essentially eliminate the possibility of any quark 
nuggets surviving till the present epoch. However, Madsen et al. [UHl then point out that 
few important effects can considerably reduce the original estimation given in ref . |H5j . 
The first important effect is related to the deficiency of u and d quarks (in contrast 
with s quark) in the surface area. This leads to the suppression of the evaporation 
rate such that B > 10 46 can be stable against evaporation In this calculation the 
penetrability of the phase boundary was assumed to be near 100%. This assumption 
was questioned in [HJJ and [SB] where it was demonstrated that nuggets with B > 10 43 
[HTJ(-B > 10 39 according to ref. jsJBj) could survive the evaporation even if the first effect 
( described above and which led to 10 -6 suppression, see [SB] for details) is neglected. 
As discussed in [22]- [HE] the limit on B may be further reduced by reabsorption. All 
these effects taken together suggest that nuggets with B > 10 30 are not ruled out and 
can survive the evaporation ^Hj- 

Our original remark here is: along with the suppression effects discuseed above, we 
have two additional effects which may further reduce the evaporation rate. 

Indeed, the core of the axion domain wall as discussed inJT],P2] has a QCD sub- 
structure with a typical scale > lGeV. It is quite obvious that this sub-structure 
certainly reduces the penetrability of particles from inside to outside, and therefore, it 
suppresses the evaporation rate. Also, the baryon charge in superconducting phase is 
in the form of the diquark condensate rather than in form of free quarks discussed in 
the previous analysis [33] -[HE]- This fact may also considerably reduce the evaporation 
rate because it requires the breaking of the Cooper pair before the evaporation becomes 
possible. This effect certinly increases the effective binding energy and decreases the 
evaporation rate. It is difficult to make a precise estimate of these effects at the moment, 
due to the many compications discussed earlier^] as well as many additional difficulties 
mentioned above. However, we believe, it is fair to say that the QCD balls with B > 10 32 
as discussed in the previous section, can safely survive the evaporation, and therefore, 
the possibility seems worth exploring. 

Now, we wish to estimate the absolute value for the dark matter number density 
noM assuming that the nonbaryonic dark matter is actually the QCD balls. In this 
case, at the QCD phase transition at T ~ T c soon after the QCD balls are formed, udm 
can be estimated as follows, 

q27T 2 q 777. AT . . 

n DM ~ 5 • 10- 9 — gXj^i (21) 
where we used the known magnitudes for the baryon to photon ratio, %/n 7 ~ 5 • 10 -10 , 
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and the dark matter to baryon ratio, Qdm/^b — 10- Numerically, for the baryon charge 
B ~ 10 32 and effective massless degrees of freedom, ~ 10 the estimation (j2Tj) leads 
to 

rT c = n D \^T c ~ 3.5 ■ 10 13 , r ~ 10 cm, (22) 

where r has the physical meaning of an average distance between QCD-balls after they 
formed. As expected, average distance r is much smaller than the horizon radius R^ 00 
at the QCD phase transition, r ~ 10~ 5 R^ CD . It is quite remarkable that r is much 
larger than the size of the QCD-ball, see eq. (|2*|) . such that QCD-balls become well 
separated soon after they formed. Besides that we expect that the QCD ball size should 
be related, through dynamics, to the correlation length £ ~ m~ l of the original axion 
field. We also expect that the spatial extend of a typical closed wall at the instant of 
formation has the same order of magnitude £ [23 12H] ■ Initial size of a closed wall ~ £ 
eventually (after some shrinking as a result of tension, and after some expansion as a 
result of evolution of the Universe) determines the size of the QCD-balls. However, the 
dynamics of this transition is quite complicated, and we are not able to derive a relation 
between initial domain wall size distribution and QCD-ball size distribution at the later 
stage. Close numerical values for the QCD ball size and £ ~ m" 1 also suggest that these 
parameters are related somehow. Therefore, it is at least possible, that the decay of the 
axion domain wall network may result in formation of the QCD-balls with their nice 
properties discussed in this work. 

To conclude: we advocate the idea that the QCD-balls could be a viable cold dark 
matter candidate which is formed from the ordinary quarks during the QCD phase 
transition when the axion domain walls form. As we argued the system in the bulk may 
reach the critical density when it undergoes a phase transition to a color superconducting 
phase in which case the new state of matter representing the diquark condensate with 
a large baryon number B becomes a stable soliton-like configuration. The scenario is 
no doubt lead to important consequences for cosmology and astrophysics, which are 
not explored yet. In particular, some unexplained events, such as Centauro events, or 
even the Tunguska-like events (when no fragments or chemical traces have ever been 
recovered), can be related to the very dense QCD balls. The recent detection |2U] of 
two seismic events with epilinear (in contrast with a typical epicentral ) sources may 
also be related to the very dense QCD balls. Also, the "missing" baryons in Galaxy 
Clusters [39 may also be related to the QCD balls. Finally, the cuspy halo problem in 
dwarf galaxies might be related to the unstable cold dark matter |40j . which, again, 
could be related to the QCD balls discussed in this work. Indeed, as we mentioned in 
the Introduction, if the QCD ball size exceeds the critical value, it becomes metastable 
(rather than stable) configuration. The life time of these metastable QCD balls could 
be very large. Therefore, they could serve as decaying dark matter particles suggested 
injUj. 

Therefore, the "exotic", dense color superconducting phase in QCD, might be much 
more common state of matter in the Universe than the "normal" hadronic phase we 
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know. More than that: one can present some arguments |41j to support the idea that the 
observed in nature asymmetry between baryons and antibaryons may also be originated 
from the same physics during the QCD phase transition. In this case the antimatter is 
hidden inside of the anti- QCD balls in the form of the diquark condensate similar to 
the QCD ball case. One could naively think that such a scenario is in contradiction with 
observations on absence of antimatter around us. However, such a conclusion would be 
very premature one due to the specific interaction features of the matter in hadronic 
phase with the matter in color superconducting phase. Namely, if the energy of the quark 
which hits the anti-QCD ball is smaller than the superconducting gap A, the quark can 
not penetrate into the bulk, break the Cooper pair and excite the internal degrees of 
freedom of the system. Rather it will be reflected jH]. A similar property is commonly 
known as the "Andreev Reflection" in the literature on conventional superconductivity. 
Therefore, at low energies, the anti-QCD balls behave as QCD balls with respect to the 
interaction with environment, and there is no contradiction with known constraints on 
such kind of anti matter in our Universe. 

In this case, without fine tuning of parameters, one can easily understand the 
relation between Qdm ~ &>b which both originated at the same instant. As is known 
this ratio is very difficult to understand if these quantities do not have the same origin. 

In conclusion, qualitative as our arguments are, they suggest that the dark matter 
could be originated at the QCD scale. 
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Appendix 

The main goal of this Appendix is to argue that the probability of the absorption of 
axion by emitted quark is extremely small. In this case our criteria of stability of the 
QCD balls, see eq.([10|h which neglects the axion domain wall energy contribution, is 
justified. Such a treatment of the problem essentially implies that we impose a weaker 
condition of metastability (rather than a stronger condition of the absolute stability) on 
the QCD balls. In different words, we assume that the energy can not be transformed 
from the axion domain wall to the quark which is about to leave the system. In what 
follows we make some estimates which support this assumption. Indeed, as we shall see 
in a moment this probability is exceedingly small due to the very small nucleon-axion 
coupling constant, ~ (1(T 13 - 1(T 15 ). 

We start from the estimation of the probability to absorb the axion from the axion 
domain wall background field a(z) by an elementary excitation \ip in > with mass A 
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proportional to the gap. This elementary excitation carries the unit baryon chage in 
superconducting phase. We assume that the final state is represented by the wave 
function of constituent quark < ip ou t\ with mass ~ mjy/3 ( hadronic phase). We take a 
simpe expression for the axion-quark interaction to be y 3 -ipl ut a(z)ipin, where m q is the 
current mass quark of order few MeV, and ~ 1 is the axion domain wall background 

Ja 

field. The precise expression for the domain wall profile function a(z) is known, however, 
in our estimate we shall use a simple expression a(z) ~ f a e~ maZ i n order to emphasize 
that the magnitude of the axion field vanishes at infinity and the typical scale where 
axion field varies is m~ l . We also assume that the quark has a trajectory z = vt with 
velocity v close to the speed of light. In this case the time-dependent interaction takes 
the form ~ m q e~ mat , and the probability for the transition can be estimated from the 
dimensional arguments as follows, 

W ~ K r dte™-™** + h.c.\ 2 ~ 1 m o qm \ | 2 , (23) 

Jo m i + ^ 

where u is the energy difference between \ip in > and |?/w > states. It is impotant to 
note that, typically u ~ (100 — 200)MeV is large, and therefore, the probability 
is very small. We neglected many factors in estimate (|23|) . In particular, we neglected 
the momentum dependence of \ip in > and \ip ut > states; the mismatch between these 
momenta would bring an additional suppression to (J23|) . and we ignore this effect at 
the moment. It is easy to understand the source of the suppression in eq. (j2~3*)) : the 
probability for a considerable excitation ~ u of the system by a smooth field with a 
typical correlation scale m~ l is very small. 

In order to derive a total number of events of absorption Wtot one should multiply 
the expression (|2 3)1 by an additional factor describing a total number of elementary 
quark excitations close to the surface of the system such that they can leave the system 
without re-scattering. This requirement (to be close to the surface of the QCD ball) is 
important because the distance from the surface should not exceed the mean free path. 
Otherwise, the quark even if it absorbs the axion, would not be able to leave the system. 
Assuming the thermodynamical equilibrium at temperature T soon after the formation 
of the QCD balls, we can estimate this factor as follows ^-g*T 3 exp(— ^)47ri? 2 ^, where 
£ is the mean free path which we estimate to be 1/T. Our final expression for the total 
probability of absorption of the axion (while the temperature is of order T) is estimated 
to be 

. m a m a l9 27r 2 o 9 , A x . , 

Wtot ~ Hrt^ 7T 9*TH<kR 2 vc? -- , 24 
m A a + uj z 45 1 

where we neglected many additional suppression factors, such as factor 1/6 describing 
the probability for the quark to move in the direction pointing off the center of the QCD 
ball. Numerically, even if we neglect the factor exp(— =) in eq. ()24j) . the probability 
is already quite small, W to t < 10 -3 for the typical values of g* ~ 10 , T ~ 0.6A. 
When temperature becomes considerably smaller than T, the probability of absorption 
diminishes due to the small number of excitations, ~ exp(— ^). This late epoch of 
evolution can be ignored. Also, one should keep in mind that the quark excitations 
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are not supported in the hadronic phase due to the confinement. Therefore, one should 
have three quarks (or quark and diquark pair from the condensate) to be organized 
in a color singlet state such that it can propagate in the hadronic phase. It definitely 
gives an additional suppression which we even did not try to estimate: the suppression 
factor (f2"4"j) is already sufficiently strong for our purposes. Therefore, our treatment of 
the problem when we use a weaker condition of metastability, see eq. (II Oj) . rather than 
a stronger condition of the absolute stability, is justified. 
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